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An application of the new formulation of the eigenchannel method [R. Szmytkowski, Ann. Phys. 
(N.Y.) 311, 503 (2004)] to quantum scattering of Dirac particles from non-local separable potentials 
is presented. Eigenchannel vectors, related directly to eigenchannels, are defined as eigenvectors of 
a certain weighted eigenvalue problem. Moreover, negative cotangents of eigenphase-shifts are in- 
troduced as eigenvalues of that spectral problem. Eigenchannel spinor as well as bispinor harmonics 
are expressed throughout the eigenchannel vectors. Finally, the expressions for the bispinor as well 
as matrix scattering amplitudes and total cross section are derived in terms of eigenchannels and 
eigenphase-shifts. An illustrative example is also provided. 
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I. INTRODUCTION 

Recently, Szmytkowski [1] proposed a new formula- 
tion of the eigenchannel method for quantum scatter- 
ing from Hermitian short-range potentials, different from 
that presented by Danos and Greiner [2]. Some ideas 
leading to this method were drawn from works on elec- 
tromagnetism theory by Garbacz [3] and Harrington and 
Mautz [4]. This method was further extended to the 
case of zero-range potentials for Schrodinger particles by 
Szmytkowski and Gruchowski [5] and then for Dirac par- 
ticles by Szmytkowski [6] (see also [7]). 

On the other hand, it is the well-known fact that sep- 
arable potentials, since they provide analytical solutions 
to the Lippmann-Schwinger equations [8] , have found ap- 
plications in many branches of physics, both in the non- 
relativistic and relativistic cases [9]. (It should be noted 
that much larger effort has been devoted to the separa- 
ble potentials in the non-relativistic regime.) Especially, 
their utility was confirmed in nuclear physics by success- 
ful use for describing nucleon-nucleon interactions [10]. 
Moreover, methods allowing one to approximate an ar- 
bitrary non-local potential by a separable one are known 
[11]. " 

In view of what has been said above, it seems inter- 
esting to pose the question: how does the new method 
apply to quantum scattering from non-local separable po- 
tentials? Partially, the answer was given by the au- 
thor by applying the method to quantum scattering of 
Schrodinger particles from separable potentials [12]. In 
the present contribution, we extend considerations from 
[12] to the case of Dirac particles. 

This paper is organized as follows. In Section 2 some 



facts and notions from the theory of potential scattering 
of Dirac particles (see [13]) are provided. In Section 3 we 
concentrate on the special class of non-local potentials, 
namely, separable potentials. In this context, expres- 
sions for the bispinor as well as matrix scattering ampli- 
tudes are provided. Section 4 contains main ideas and 
results. Here, we define eigenchannel vectors, directly re- 
lated to eigenchannels, as solutions to a certain weighted 
eigenproblem. Moreover, we introduce eigenphase-shifts, 
relating them to eigenvalues of this spectral problem. 
Within this approach, we also calculate expressions for 
the scattering amplitude and the average total cross sec- 
tion. In Section 5, scattering from a rank one delta-like 
separable potential is discussed as an illustrative exam- 
ple. The paper ends with two appendices. 



II. QUANTUM SCATTERING OF DIRAC 
PARTICLES FROM NON-LOCAL POTENTIALS 



Let us assume that a free Dirac particle of energy E 
(with \E\ > mc 2 ) described by the following monochro- 
matic plane wave 



&(r) = (r\k lXl ) = l/iO^e**-*, 



where 



Xi 



Vl + £ 2 V £(T ' k * Xi 



E — mc 2 



E 



(2.1) 



(2.2) 



(2.3) 
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is being scattered from a non-local potential given by a 
kernel V(r,r'), which in general may be a 4 x 4 matrix. 
In the above equation, Xi stands for a normalized pure 



2 



spin-i state belonging to C 2 . Orientation of the spin in 
R 3 will be denoted by Vi and is related to Xi by v% — 
x\ a Xu where cr is a vector consisting of the standard 
Pauli matrices, i.e., 



1 

1 



-i 

1 



1 

-1 



(2.4) 



Moreover, = hki is a momentum of the incident par- 
ticle and k denotes the Dirac wave number and is given 

by 



k = sga(E)] 



E 2 — [mc 2 Y 



c 2 n 2 



(2.5) 



Thereafter, we shall consider only Hermitian potentials, 
i.e., those with kernels obeying V(r, r') = V^(r', r). 

For this scattering process we may write the 
Lippmann-Schwinger equation [8] of the form 

V>(r) = cpi(r) 

-/ dV f rf 3 r"G(S,r,r')V(r',r")^(r"). (2.6) 

JR 3 JR 3 

Function G(E, r, r') appearing above is the relativistic 
free-particle outgoing Green function given by 



1 



G(E, r, r') = (-icha • V + (3mc 2 + El 4 ) 



a ik\r-r'\ 



|r - r'| 
(2.7) 

and formally is a kernel of the relativistic outgoing Green 
operator defined as 

G(E) =lim[W -E-ie}-\ (2.8) 

with Ho = —icha ■ V + /3mc 2 being a Dirac free-particle 
Hamiltonian. Here 



cr 
cr 



1 2 
-1 2 



Is 



1 
1 

,.( 2 - 9 ) 

and I4 = I2 <S> 1-2- It is worth noticing that within the 
relativistic regime the Green function (2.7) is a 4 x 4 
matrix. 

For purposes of further analysis, it is useful to intro- 
duce the following projector: 



V(k) 



cha ■ k + /3mc 2 + El 4 
2E ' 



(2.10) 



which, as one can immediately infer, may be decomposed 
in the following way 



r(k) = e+(k)e;(k) + e_(k)e t _(k) 

1 / 1 2 ecr • k 
1 + e 2 \so--k e 2 l 2 , 

with 6±(k) being defined as 

6±(k)= VTT^(-'^ ± 



(2.11) 



(2.12) 



Spinors 0± constitute an arbitrary orthonormal basis in 
C 2 , i.e., 6% = 6 st (s,t = -,+) and £+ = _ 6„6\ = 1 2 . 
What is important for further considerations, the matrix 
(2.10) possesses the obvious property that 7 , (k)9±(k) = 
0±(k) and therefore 



P(k i )I7 < (ki) = UiQm 



(2.13) 



We shall be exploiting this property in later analysis. 

Considering scattering processes we usually tend to 
find expressions for a scattering amplitude and various 
cross sections. To this aim we need to find an asymp- 
totic behavior of the relativistic outgoing Green function. 
From Eq. (2.7), using the projector (2.10), we have 



G(E,r,r') 



E 



■ikr 



2irc 2 H 2 



,-ik 



(2.14) 



where kf = kr/r is a wave vector of the scattered parti- 
cle. Notice that due to the fact that we deal with elastic 
processes |kj| = |k/| = k. After application of Eq. (2.14) 
to Eq. (2.6), we obtain 

V>(r) r ~°° asymp0i(r) + A fi , (2.15) 

r^oo f 

where Afi is the bispinor scattering amplitude and is 
defined through the relation 



E 



2nc 2 h 2 



V{k f ) 



f d 3 r' f d 3 r"e-^ r 'V(r',r")?A(r") (2.16) 

JR 3 JR 3 



and, in general, is of the form 



err ■ k fXf 



(2.17) 



Here x/ is a spinor transformed from the initial spinor \i 
by the scattering process. Vector Uf — (X/ cr X/)/(X/X/) 
responds for an orientation of the spin of the scattered 
particle. Therefore let us assume that there exist a ma- 
trix, such that Xf — ^fiXi- Then it is easy to verify that 
the bispinor scattering amplitude may be written in the 
form 

A fi = AfiUiiki), (2.18) 
where the matrix Afi is related to s/fi by 



■Afi 



fi 



sjafficr ■ ki 



1 + e 2 \ ecr • kfs/fi e cr ■ kfg/fi(T ■ ki 



(2.19) 

Henceforth matrices Afi and srffi will be called the ma- 
trix scattering amplitudes. The differential cross section 
for scattering from the direction kj and the spin arrange- 
ment v\ onto kf and i>»f is defined as 



Aa_ 

dsTT 



XfXf = xWfi^fiXi 



(2.20) 
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Subsequently, after integration the above over all the di- 
rections of k f , we arrive at the total cross section 



o-(ki,i/i) = i> d 2 k f Xf Xf- 



(2.21) 



Finally, averaging over all directions of incidence kj and 
the initial spin orientation C>i 1 one finds the average total 
cross section 



<Jt(E) 



1 



. ., d% d> d 2 Ui J> d 2 k f X )xf- (2-22) 

Jin Jiir Jin 



Obviously all the mentioned cross sections may be ex- 
pressed in terms of all the scattering amplitudes Afi, 
Afi, and sifi. 



III. SPECIAL CLASS OF NON-LOCAL 
SEPARABLE POTENTIALS 

In this section we employ the above considerations to 
the special class of non-local separable potentials. As 
previously mentioned, such a class of potentials allows to 
find solutions to the Lippmann Schwingcr equations in 
an analytical way. 

Consider the following class of potential kernels: 



V(ry) = 5> M iv(r)ut(r') 



(3.1) 



where it is assumed that in general /i may denote the 
arbitrary finite set of indices, i.e., /i — {hi, ■ ■ ■ , Mfc} an d 
all the coefficients lo^ different from zero. Functions u M (r) 
are assumed to be four-element columns. 

Substitution of Eq. (3.1) to Eq. (2.6) leads us to the 
Lippmann-Schwinger equation for the separable poten- 
tials: 

V(r) - &(r) - V w„ / dV G(E, r, r')u^(r') 

x / d 3 r" U t(r>(r"), (3.2) 

which may be equivalently rewritten as a set of linear 
algebraic equations. Indeed, using the Dirac notation 
one finds 

[<W + (u.|G(£)KK] = (u v |0i). (3.3) 

For further convenience we introduce the following nota- 
tions 

/ (ui\ip) 
(u\<p) = <u 2 |^) 

<^|u> = <u|^> + = (<^| Ul > <^|u 2 > ...)• (3.4) 



Consequently, we may rewrite Eq. (3.3) as a matrix equa- 
tion (1 + GQ)(u|f/') = (u|0i) or equivalently as 



(u|^) = (l + GQ)- 1 (u|^), 



(3.5) 



with G being a matrix composed of the elements 
(u ! y|G(i?)|u A1 ) and ft = diagju^]. Similarly, substitut- 
ing Eq. (3.1) to Eq. (2.16) and again using Eq. (2.10), 
we arrive at the bispinor scattering amplitude for the 
separable potentials in the form 



-E 



x f dV ut (r')V(r'), (3.6) 

which, by virtue of Eqs. (2.11) and (3.5), reduces to 
-E + 

'' ^e s (k / )(k / s |u)ft(l + Gft)" 1 (u|^) 



27TC 2 ft 2 



(3.7) 

and, utilizing the fact that for all invertible matrices X 
and Y the relation (XY) _1 = Y _1 X _1 is satisfied, finally 
to 

Afi = 2^ ^ Qs(k f )(k f 8 s \u) (ft- 1 + G)- 1 <u|^>. 

(3.8) 

Subsequently, using the fact that (2.13), we obtain the 
bispinor scattering amplitude in the following form 



A fi 



-E 



2irc 2 h 2 



J2 e s (k f ){k f e s \u) 



s,t= 



x(ft- 1 + G) 1 (u|k i e t )eJ(k i )i; < (k i ), (3.9) 



which, after comparison with Eq. (2.18), gives the for- 
mulae for the 4x4 matrix scattering amplitude: 



Aft 



-E 



2irc 2 H 2 



J2 e s (k / )(k / ^|u) 



x (ft-^G) 1 (u|k i ^ t )e t t (k i ), 



(3.10) 



and finally, after straightforward movements, for the 2x2 
matrix scattering amplitude as 



-E 



2nc 2 h 2 



J2 e s (M s |u)(ft- 1 + G) '(ulMt)^- 

s,t=- 

(3.11) 



IV. THE EIGENCHANNEL METHOD 

Now we are in position to apply the eigenchannel 
method proposed recently by Szmytkowski [1] to scat- 
tering of the Dirac particles from potentials of the form 



4 



(3.1). As we shall see below, such a class of potentials al- 
lows us to formulate this method in a simplified algebraic 
form. 

We start from the decomposition of the matrix ft -1 + 
Gd into its Hcrmitian and non-Hcrmitian parts, i.e., 



ST 



A -MB, 



(4.1) 



where matrices A and B are defined through relations 

A = ft- 1 + i(G + Gt), B =1(G-Gt). (4.2) 

It is evident from these definitions that both matrices A 
and B are Hermitian. Moreover, utilizing the fact that 



ike ik \ r 



■ik\i 



(4.3) 



where g = r — r', the straightforward calculations lead 
us to their matrix elements of the form 



(r) 



ichka ■ JL yi (k\Q\) + (/3mc 2 + E)y (k\g\) 

\q\ 



(4.4) 



and 



JR 3 JR3 



d 3 r'ut(r) 



ichka. ■ JL^{k\Q\) + {I3mc 2 + E)j (k\g\) 
\Q\ 



(4.5) 



where jo (z), ji{z), yo{z) and yi(z) are, respectively, the 
Bessel and Neumann spherical functions [14]. Recall that 
in general jo(z) = (sinz)/z, yo(z) = (cos z)/z and 



, cosz sin z 
3i{z) = + 



sin z 



cos z 



■ (4-6) 



The main idea of the present paper, adopted from [1] , is 
to construct the following weighted spectral problem: 



AX 7 (£) = X 7 (E)BX 1 (E), 



(4.7) 



where X 1 (E) and X 7 (E) is, respectively, an eigenvector 
and an eigenvalue. Thereafter the eigenvectors {X 1 (E)} 
will be called eigenchannel vectors. They are directly 
related to eigenchannels defined in [1] as state vectors. In 
fact, they constitute a projection of eigenchannels onto 
subspacc spanned by u p (r). 

Using the fact that matrices A and B are Hcrmitian 
and, as it is proven in Appendix A, the matrix B is posi- 
tive semi-definite, one finds that the eigenvalues {A 7 (i?)} 
are real, i.e., X*(E) — A 7 (_E). Moreover, the eigenchan- 
nels associated with different eigenvalues obey the or- 
thogonality relation 

xUe)BX,(E)=0 (X y (E)^X 1 (E)). (4.8) 



In case of degeneration of some eigenvalues one may al- 
ways choose the corresponding eigenvectors to be orthog- 
onal according to the above relation. Then, imposing 
the normalization XUE)BX^(E) = 1, one obtains the 
following orthonormality relation 



xUe)bx^(e) 



(4.9) 



From Eqs. (4.7) and (4.9) one infers that the eigenvalues 
{X 1 (E)} may be related to the matrix A as follows 



Ky(E) = XhE)AX^(E). 



(4.10) 



Similar reasoning may be carried out employing the ma- 
trices A and ft -1 + G. Indeed, after algebraic manipula- 
tions we arrive at 

x\,{E)kX 1 {E) = X 1 {E)8 1 , 1 , 
A^Xft- 1 + G)X 1 {E) = [i + A 7 (£)]<V 7 , (4.11) 

and A 7 (£) = X^E)^ 1 + G)X 1 {E) - i. Since the 
eigenchannels {A 7 (_E)} are the solutions of the Hermi- 
tian eigenvalue problem, they may satisfy the following 
closure relations 

Y J X,{E)X\{E)B = t, 

7 

Y,^\E)X 1 {E)X\{E)k=l, (4.12) 



and 



A 7 (£) 



X^{E)XhE)(Sr x + G) = 1, (4.13) 



where 1 is an identity matrix, which dimension depends 
on the dimension of the matrix G. For purposes of further 
analyzes the above closure relations are assumed to hold. 
Below, we employ the above reasoning to the derivation 
of the scattering amplitudes. From Eq. (4.13) one de- 
duces that 



(ft-i+cr^E 



i 



i + A 7 (£) 



X,(E)Xl(E). (4.14) 



After substitution of Eq. (4.14) to Eq. (3.10) and rear- 
ranging terms, we have 



—E 



i 



2Trc 2 h 2 ^i + X 7 (E) 



J2Qs(k f ){k f e s \u)x y (E) 



^x^EX^e^eli^). 



(4.15) 



Let us define the following angular functions 
3^ 7 (k) 



V^^£© S (k)(k^|u)X 7 (£), (4.16) 



■5 



hereafter termed the eigenchannel bispinor harmonics. 
The functions {3^ 7 (k)} are orthonormal on the unit 
sphere (for proof, see Appendix B), i.e., 

I d 2 kyl,(k)y 7 (k) = s Yl . (4.17) 

J4-7T 

Application of Eq. (4.16) to Eq. (4.15) yields 



,M 7 (S) 



k 



sin S^y^yUki), (4.18) 



where {8 7 (E)} are called eigenphase- shifts and are re- 
lated to {X 1 (E)} according to 



X 7 (E) = - cot 5 7 (E). 



(4.19) 



Similar considerations may be carried out for the 2x2 
matrix scattering amplitude &/fi. Indeed, in virtue of Eq. 
(2.19) we may rewrite Eq. (3.11) in the form 

*>« = yE ea,(E)sin M £ ) T T(k/)TiW, (4-20) 

7 

where the angular functions {T 7 (k)}, hereafter called 
eigenchannel spinor harmonics, are defined as follows 



r 7 (k) = 



Ek 



87T 2 C 2 ft 2 



^2e a (k0 B \u)X^(E). (4.21) 



Moreover, they are orthogonal on the unit sphere (for 
proof, see Appendix B) 



d 2 kTL,(k)T 7 (k) = S, 



7 ' 7 , 



(4.22) 



4tt 



and, as one can verify, are related to the eigenchannel 
bispinor harmonics {3^ 7 (k)} via the relation 



3> 7 (k) 



1 



T 7 (k) 



yTTT 2 " V ecr ' kT 7( k ) 



(4.23) 



Now we are in position to compute scattering cross- 
sections. Substitution of Eq. (4.20) to Eq. (2.20) and 
integration over all directions of scattering k/ , by virtue 
of relation (4.22), yields 



cr(kj,i/j) 



167T 2 

k 2 



J2sm 2 6,(E) xlT 7 (k t ) . (4.24) 



To compute the total cross-section averaged over all ar- 
rangements of spin of the incident particle, we have to 
notice that the projector onto the pure state \i mav 
be written as XiXi = (l/ 2 )[l2 + v>i ■ cr] with = 1. 
Therefore, substituting of the above to Eq. (4.24) and 
averaging over all directions of i/j , we arrive at 

^) = ^E sin2 ^)n( k *) T 7(k*)- (4.25) 



Finally, averaging the above scattering cross-section over 
all directions of incidence k^, again by virtue of Eq. 
(4.22), we get the total cross-section in the form 



2vr 



(4.26) 



It should be emphasized that all the above considerations 
respecting scattering cross-sections may be repeated us- 
ing the eigenchannel bispinor harmonics {^(k)} instead 
of the eigenchannel spinor harmonics {T 7 (k)}. The sig- 
nificant difference is that then the integrals over ky and 
ki need to be calculated using relation (4.17) instead of 
(4.22). 



V. EXAMPLE 

We conclude our considerations providing here an illus- 
trative example concerning the scattering from a spher- 
ical shell of radius i?, centered at the origin of the co- 
ordinate system. Due to the assumption of non-locality 
of potentials under consideration, we shall simulate this 
process by using a potential of the form 



V(r,r') = wv(r)v(r')l 4 , v(r) = 



1 S(r - R) 



(5.1) 

where w^O. Notice that the potential defined above is 
the special case of that proposed recently by de Prunele 
[15] (see also [16]). Scattering of the Dirac particles from 
delta-like potentials was also studied e.g. in Refs. [17, 
18] . However, in these papers the authors considered only 
local potentials and not non-local ones. 

At the very beginning, we need to bring the potential 
(5.1) to the previously postulated form (3.1). To this 
aim, let ei and e2 constitute a standard basis in C 2 , i.e., 
ei = (1 0) T and e 2 = (0 1) T . Moreover, let e y - = e 2 <£> ej 
and then by virtue of the fact that I4 = J2ij=i e ij e tj> 
we may rewrite (5.1) as 

2 

V(r,r') = u> J2 ( r )> M r ) = «( r )ey- 

(5.2) 

Now, we are in position to compute the matrix G. Using 
Eqs. (2.7) and (5.1), after straightforward integrations 
we have 



G = iki {kR)h [ +\kR) 



J?+l 2 








(5.3) 



where n± = (E±mc 2 )/c 2 h 2 and h^ + \z) = j a (z) + iy (z) 
is the spherical Hankel function of the first kind. Hence, 
by the definitions given in Eq. (4.2), we find that the 
explicit forms of matrices A and B are 

o- 1 - kj a {kR)y {kR)n + }l 2 \ 
[uj- 1 - kj (kR)y (kR)r]-}l2 ) 

(5.4) 
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and 



B = kf (kR)(^ 







(5.5) 



According to the method formulated in Sec. IV, we may 
construct the following spectral problem 

AX y (E) = X 1 (E)BX y (E) (7=1,2,3,4), (5.6) 

which, as one can easily verify, has two different eigen- 
values 



\±(E) = 



w 1 - kj (kR)y (kR)i] d 



(5.7) 



and respective eigenvectors 
\E) = 

- > [2)) {E) 



kr] + jo(kR) 
1 



ei 



e 2 



ei (2 ), 
>ei( 2 ). 



/kr)-j (kR) 

Then, using Eq. (4.16) and by virtue of the fact that 



(5.8) 



(kx|u) 



47T 



Jo(kR) 



1 + e 2 



(5.9) 



we arrive at the four cigenchannel bispinor harmonics 
{y^ (k)} in the form 



3><. 1(2 »(k) 



2l(2) 



and 



5,(1(2)) (k) 



yjAir(l + e 2 ) V £(T ' ke i(2) 



err • ke 1(2) 
ei( 2 ) 



V4tt(1 + e 2 ) 



(5.10) 



(5.11) 



Then, by virtue of Eq. (4.21), one obtains the eigenchan- 
nel spinor harmonics {T 7 (k)} in the form 



-(1(2)) 



(k) 



ei (2 ), 



-(1(2)) 



(k) 



1 



a ■ ke 1(2) . 



/47T V47T 

(5.12) 

The latter may be equivalcntly obtained combining Eqs. 
(4.23) and (5.12). Moreover, as one may easily verify, 
functions given by Eqs. (5.10) and (5.12) are orthonor- 
mal, respectively, in the sense (4.17) and (4.22). 

Before we find an expression for total cross section, we 
compute the scattering amplitude. Since, as shown in 
Sec. II, the bispinor and both matrix scattering ampli- 
tudes are mutually related, we restrict our considerations 
to the 2x2 scattering amplitude. Thus, combining Eqs. 
(4.20), (5.5), and (5.12) we obtain 



s/ fi = -jlikR) 



Jk JO (kR)h ( + \kR) + (u JV+ )-i 
(<r-k / )(<r-k i ) 



ikj (kR)h ( +) (kR) + (wj?-)" 



(5.13) 



Finally, substitution of Eqs. (5.8) and (5.12) to Eq. 
(4.24) with the aid of Eq. (4.19) yields 

47T 4 

(j{ki,Vi) = —j {kR) 



l(kLor, + )-i - j Q (kR)y a (kR)] 2 + ^(kR) 



+ 



[(fc^.)- 1 - 3o (kR)y (kR)} 2 +J$(kR) 



}■ 



(5.14) 



Here it is evident that cx(kj, i/,) = <r(kj) = a t (E). 

In order to illustrate the obtained results, the eigen- 
phaseshifts for two different values of u), derived from 
Eqs. (4.19) and (5.7), are plotted in Fig. 1 and 2. Fig- 
ures 3 and 4 present partial a±(E) as well as total <r t (E) 
cross sections. 

It seems interesting to investigate the behavior of both 
eigenvalues \±(E) in the non-relativistic limit, i.e., when 
c — > co. From (2.5) one concludes that 



c^oo 2m 
V+ 'IP' 



n- 



o 



and therefore 



\ + (E) 



=^oo (h 2 /2m<jj) - kj (kR)y (kR) 



and 



X_(E) 



kj 2 (kR) 



sgn(w)co. 



(5.15) 



(5.16) 



(5.17) 



This means that 6 -(E) — ► nir (n e Z) in the limit of 
c — > oo. Therefore the cross section <T-(E) vanishes in 
the non-relativistic limit and in this sense it has a purely 
relativistic character leading to the fact that the reso- 
nance appearing in Fig. 4 at about 1.25mc 2 is purely 
relativistic effect. 

One sees that in the non-relativistic limit the cross 
section (5.14) reduces to 



crt(E) 



47T 

k 2 



Jo(kR) 



[(h 2 /2mkuj) - j (kR)y (kR)] 2 + j*(kR) 



.(5.18) 



The above cross section may also be obtained using non- 
relativistic formulation of the present method given in 
Rcf. [12]. 



VI. CONCLUSIONS 

In this work, an application of the recently proposed 
eigenchannel method [1] to the scattering of Dirac par- 
ticles from non-local separable potentials has been pre- 
sented. Application of such a particular case of the non- 
local potentials reduces naturally the general weighted 
eigenvalue problem stated in Ref. [1] to its matrix coun- 
terpart given by Eq. (4.7) leading to the definition of 
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FIG. 1: Behavior of eigenphaseshifts 8+(E) (solid curve) and 
6 -(E) (dashed curve) as functions of energy E (in units of 
mc 2 ) for ui = — h 3 /m 2 c and R = h/rnc. The eigenphaseshift 
S+(E) has been constrained to the range [— it/2, n/2]. 
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FIG. 2: Behavior of eigenphaseshifts 5+(E) (solid curve) and 
8 -(E) (dashed curve) as functions of energy E (in units of 
mc 2 ) for lu — —5h 3 /m 2 c and R = h/mc. Both eigenphase- 
shifts have been constrained to the range [—n/2, n/2]. 



eigenchannel vectors. Using the notion of the eigenchan- 
nel vectors the definitions of eigenchannel spinor as well 
as bispinor harmonics have been given. The latter pro- 
vide us with the formulas for scattering amplitudes simi- 
lar to that well-known for central potentials generalizing 
them at the same time to the case of non-local separable 
potentials. 

The general considerations have been extended with 
an illustrative example in which the Dirac particles are 
scattered from non-local, delta-like potential. In this par- 
ticular case, the general eigenvalue problem (4.7) become 
just a 4 x 4 matrix equation and therefore is easily solv- 
able (notice that in the case of non-relativistic scatter- 
ing it would be just a one-dimensional problem). The 
eigenvalues of this problem are two-fold degenerated and 
therefore give two different eigenphasc-shifts from which 
one has a purely relativistic character in the sense that it 
tends to nn (n € Z) whenever c — > oo giving no contribu- 
tion to total cross sections in non-relativistic limit. One 
sees also that even such a simple example of non-local 
potentials may lead to some resonances (see Fig. 4). 

The next step in our considerations will be to inves- 
tigate the applicability of the new formulation of the 
eigenchannel method in the case of inelastic scattering 
from separable potentials. Moreover it seems also inter- 
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FIG. 3: Partial <J+(E) (dashed curve), <J-(E) (dotted curve), 
and total <Jt(E) (solid curve) cross sections (all in units of R 2 ) 
as functions of energy E (in units of mc 2 ) for u> = — ft 3 /m 2 c 
and R — h/mc. 
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FIG. 4: Partial a+(E) (dashed curve), <J~(E) (dotted curve), 
and total a t (E) (solid curve) cross sections (all in units of R 2 ) 
as functions of energy E (in units of mc 2 ) for uj — —5h 3 /m 2 c 
and R — h/mc. 



esting to investigate the applicability of the method to 
the other, more complicated examples of separable po- 
tentials. 
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APPENDIX A: POSITIVE SEMIDEFINITENESS 
OF THE MATRIX B 

The proof follows the suggestions of Szmytkowski [19]. 
Positive scmidcfmiteness of the matrix B means that the 
inequality 

XUE)BX 7 (E) > (A.l) 
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is satisfied. To prove the above statement let us notice 
that 

I d 2 ke^ e {cha-k + (3mc 2 + El 4 ) = 4tt 



icKkji{k\o\)a ■ f-r + (fimc 2 + El 4 )j (k\g\) 



(A.2) 



where g = r r'. Then using Eq. (2.10), we may rewrite 
Eq. (4.5) in the form 



8TT 2 C 2 fl 2 



I d 2 k [ d 3 re tk r u, t ,(r)'P(k) 

JilT JR. 3 

[ d 3 r' e-^'u^r'), (A.3) 

JR 3 



which after application to Eq. (A.l) yields 
X\(E)BX 7 (E)= " ,k 



8n 2 h 2 



<f d 2 k Y, X U E ) I d 3 re ikr ut(r)P(k) 

Jin „ Jr 3 



>0, 
(A.4) 



finishing obviously the proof. Here X 1V (E) denotes the 
vth element of the eigenchannel vector X 1 (E) and | |f2| | = 



APPENDIX B: ORHTONORMALITY OF THE 
ANGULAR FUNCTIONS y 7 (k) AND T 7 (k) 

We begin with proof for the functions y 7 (k) . Applica- 
tion of Eq. (4.16) to Eq. (4.17) with the aid of Eq. (2.11) 
and the fact that V{k) is a projector, we can deduce that 

£rf 2 k^ (km(k) = 



/ d 3 r f dVu£(r) 

Jr 3 Jr 3 

x I d 2 ke 4k - (r - r ' ) 7'(k)u AJ (r') 

Ji-K 



X^(E), (B.l) 



Comparison with Eq. (A.3) shows that the square brack- 
ets in the above contain the expression proportional to 
certain element of the matrix B. Therefore, we may 
rewrite Eq. (B.l) as 

<f d 2 ky f y (k)y 7 (k) = x f y (E)Bx 7 (E). (b.2) 

J 47T 

Finally, substitution of Eq. (4.9) to Eq. (B.2) leads us 
directly to Eq. (4.17), finishing the proof. To prove the 
orthonormality relation for the functions {T 7 (k)}, it suf- 
fices to combine Eq. (4.23) with Eq. (B.2). 
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